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Abstract. The role of magnetic fields in gamma-ray burst (GRB) flows remains controversial. The study of the early afterglow 
phases and, in particular, of the reverse shock dynamics and associated emission offers a promising probe of the magnetization 
of the ejecta. In this paper, we derive the conditions for the existence of a reverse shock in arbitrarily magnetized ejecta that 
decelerate and interact with the circumburst medium. Both constant and wind-like density profiles are considered. We show, 
in contrast to previous estimates, that ejecta with magnetization ctq > 1 are not crossed by a reverse shock for a large fraction 
of the parameter space relevant to GRB flows. Allowing for shell spreading, there is always a relativistic or mildly relativistic 
reverse shock forming in ctq < 0.3 ejecta. From this, we conclude that the paucity of optical flashes, believed to be a distinctive 
signature of a reverse shock, may be explained by the existence of dynamically important magnetic fields in the ejecta. 
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1. Introduction 

The role of magnetic fields in the GRB flow remains un- 
certain. Neutrino annihilation may be the dominant process 
that leads to the launching of the GRB flow by forming a 
fireball, i.e. a flow initially dominated by its thermal energy 
(Goodman 1986; Paczynski 1986). Alternatively, if they are 
powerful enough and develop the appropriate topology, mag- 
netic fields can efficiently extract the rotational energy from 
the central engine launching a Poynting-flux-dominated flow 

' (hereafter PDF; Usov 1992). 

In the fireball model the magnetic fields are dynamically 
unimportant. The GRB flow during the onset of the afterglow 

' phase is expected to be, at most, weakly magnetized. Though 
model dependent, the conversion of magnetic into kinetic en- 
ergy in a PDF is partial (Michel 1969; Li et al. 1992; Vlahakis 
& Konigl 2003; Giannios & Spruit 2006). The flow is ex- 
pected to have a significant fraction of its energy in the form 
of Poynting flux at large radii where there is substantial inter- 
action with the external medium. 

The study of the early phases of the ejecta interaction with 
the external medium can probe the strength of magnetic fields 
in flow. The early emission from the forward shock, for ex- 
ample, depends on the magnetization of the ejecta (Lyutikov 
& Blandford 2003; Zhang & Kobayashi 2005, hereafter ZK05; 
Genet et al. 2006). Modeling of the microphysical parameters 
of the forward and the reverse shock supports the existence of 
strong magnetic fields in the ejecta (Fan et al. 2002; Kumar 
& Panaitescu 2003; Zhang et al. 2003). Furthermore, polariza- 
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tion measurements can provide information on the structure of 
the magnetic fields (Lazzati et al. 2003; Granot & Taylor 2005; 
Mundell et al. 2007). 

The early optical afterglow lightcurves depend critically on 
the very existence of a reverse shock into strongly magnetized 
ejecta. A reverse shock may also contribute in a dominant way 
to the observed emission in later afterglow phases (Genet et al. 
2007; Uhm & Beloborodov 2007). In this paper, we investi- 
gate the parameter space of the magnetization of the flow for 
which there is a reverse shock crossing the ejecta. This parame- 
ter space is described by simple analytic expressions. We show 
that, in contrast to previous claims, if the flow contains similar 
amounts of kinetic and magnetic energy at the onset of the af- 
terglow, there is no reverse shock forming in a large fraction of 
GRBs. The observational implications of our findings are also 
discussed. 



2. The early afterglow phases 

2.1. Magnetization of the flow at the onset of the 
afterglow 

For the GRB flow to be accelerated to ultrarelativistic speeds, 
it must be launched with a high energy-to-rest-mass ratio. The 
primary form of the energy of the outflow remains unknown. 
In the fireball model (Goodman 1986; Paczynski 1986), most 
of the energy is assumed to be initially stored in thermal form 
while in PDFs (Usov 1992; Meszaros & Rees 1997b; Spruit et 
al. 2001; Van Puten & Ostriker 2001; Uzdensky & MacFadyen 
2006) in magnetic energy density. 
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Neutrino-antineutrino annihilation taking place in the po- 
lar region of a hyper-accreting central object can lead to the 
formation of a fireball (e.g., Woosley 1993; RufFert & Janka 
1998; Aloy et al. 2000; Aloy, Janka & Miiller 2005). The 
fireball is accelerated by its internal pressure gradients that 
convert its thermal energy into kinetic energy of the baryons. 
After the acceleration phase is over, the flow may consist of 
faster and slower shells that colUde with each other. These col- 
lisions lead to internal shocks that power the prompt emis- 
sion (Paczynski & Xu 1994; Rees & Meszaros 1994; Daigne 
& Mochkovitch 1998; Mimica et al. 2005; Mimica, Aloy & 
Miiller 2007). Internal shocks can lead to amplification of 
seed magnetic fields due to, e.g., the two-stream instability 
(Medvedev & Loeb 1999). The strength of these small scale 
fields and whether their survival on macroscopic scales down- 
stream of the shocks is rather uncertain (e.g. Spitkovsky 2007). 
Still, they may account for ~ 0.1 - 1% of the luminosity of the 
flow. After the phase of internal collisions, the flow expands 
and cools before it enters the afterglow phase. There, the flow 
is expected to be weakly magnetized (i.e. the magnetic energy 
density is much less than the rest mass energy density of the 
flow). 

On the other hand, magnetic fields can be eflicient in ex- 
tracting rotational energy from the central engine which may 
be an accretion disk (Blandford & Payne 1982), a rotating 
black hole (Blandford & Znajek 1977) or a millisecond mag- 
netar (Usov 1992). The resulting flow is initially Poynting-flux 
dominated. The acceleration of a PDF is more model depen- 
dent than that of the fireball model. It depends on the mag- 
netic field geometry and on the role of magnetic dissipation. 
In particular, magnetic dissipation can lead to efficient accel- 
eration of the flow by converting Poynting flux in to kinetic 
energy (Drenkhahn & Spruit 2002; Giannios & Spruit 2006). 
Magnetic dissipation can also efficiently power the prompt 
emission (Thompson 1994; Spruit et al. 2001; Lyutikov & 
Blandford 2003; Giannios 2006b; Thompson 2006; Giannios 
& Spruit 2007). Despite the model dependent aspects of MHD 
jet acceleration, a common outcome of different studies is that 
the acceleration process is not expected to be 100% efficient in 
converting the Poynting flux into kinetic flux. As a result, the 
flow remains strongly magnetized after the acceleration phase 
is over. At large radii, magnetic fields carry energy compara- 
ble to that of the kinetic energy of the baryons (Li et al. 1992); 
Drenkhahn & Spruit 2002; Vlahakis & Konigl 2003, Giannios 
& Spruit 2006) or even dominate over the kinetic energy of 
the flow (e.g. Lyutikov & Blandford 2003; Thompson 2006; 
Lyutikov 2006). 

It has to be noted that the degree of magnetization of the 
ejecta may critically depend on the likelihood of baryon mixing 
in the course of the generation of a PDF, which can only be re- 
liably determined by extremely well resolved MHD numerical 
simulations. A substantial baryon entrainment could result in a 
much smaller magnetization than predicted by the above men- 
tioned models. Furthermore, the magnetization of PDFs could 
be substantially decreased by the pair-loading resulting from 
vv-annihilation close to the central engine (see, e.g., Levinson 
& Eichler 1993; Aloy & Obergauhnger 2007). 



2.2. Initial interaction of the ejecta with the external 
medium 

If far enough away from the central engine, there is substantial 
interaction of the relativistic ejecta with the external medium. 
This interaction is beheved to result to the afterglow emis- 
sion. A characteristic difference between fireballs and PDFs 
lies in the magnetization of the ejecta at the onset of the after- 
glow phase. In fireball models, the energy of the flow is domi- 
nated by the kinetic energy of baryons. If the flow is launched 
Poynting-flux dominated, it is expected to maintain a large frac- 
tion of its energy in the form of magnetic energy at large dis- 
tance. 

Since the initial phases of the ejecta-external medium inter- 
action depend on the magnetization of the ejecta (e.g., Lyutikov 
2006), it is possible distinguish among fireball and PDF mod- 
els using early afterglow observations. A particularly promis- 
ing probe of the magnetic content of the flow comes from the 
study of the reverse shock emission (e.g. Kumar & Panaitescu 
2003; Zhang et al. 2003). In this work, we focus on the condi- 
tions for the existence of a reverse shock in arbitrarily magne- 
tized ejecta. 

We focus on the interaction of the GRB ejecta with exter- 
nal medium of constant number density n^. Similar analysis 
can, however, be appUed to different density profiles for the ex- 
ternal medium. The case of wind-like density profile n ~ r^^ 
is also discussed in Sect. 3.2. For typical GRB parameters, this 
interaction becomes substantial at distances r ~ lO'^ - lO'^ cm. 
We assume that the acceleration and collimation of the flow has 
taken place much closer to the central engine as well as the in- 
ternal dissipation mechanisms (e.g. internal shocks, magnetic 
dissipation) that are believed to power the prompt emission. 
After internal dissipation is over, the flow expands and cools 
down. Since we are interested in the afterglow phase, we con- 
sider a, radially moving, cold and relativistic shell with width 
Ao and corresponding Lorentz factor 70 » 1. 

In a radially moving magnetized flow, the induction equa- 
tion results in a component of the magnetic field perpendicular 
to the motion B± oc 1/r that drops much slower with radius 
than the parallel one: B|| oc l/r^. Since the flow expands for 
several order of magnitude from the launching to the afterglow 
region, we assume that the magnetic field is dominated by its 
perpendicular component. The magnetic content of the shell 
is conveniently parameterized by its magnetization ctq which 
stands for the Poynting-to-kinetic flux ratio at some initial ra- 
dius fQ. For an ultrarelativistic flow 



Ek 47TyopQC^ ' 
where Bq and po are the initial magnetic field strength and den- 
sity of the shell; both measured in the lab frame. Note that, 
since the shell is assumed initially cold, we neglect the contri- 
bution of the internal energy to the kinetic flux of the shell. 
The total energy of the shell is 

E = 4nRlAoiyoPoc^ + Bl/4n) = E^il + (Tq), (2) 

where the first term stands for the kinetic energy and the second 
for the magnetic energy of the sheU. In other words, the mag- 
netization (To parameterizes the fraction of the total energy of 
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the flow that is in kinetic and in magnetic form. These fractions 
are 1/(1 + <To) and (To/(1 + ctq) respectively. 

Once E, Aq, yo and cro are fixed the properties of the ejecta 
are defined. The reference values that we adopt for these pa- 
rameters are summarized in the following. In our spherically 
symmetric model E is roughly equal to the isotropic equiva- 
lent energy Ei^o of the burst and, thus E ~ Ei^o ~ 10^"'£'53 erg 
(with spread of a couple of orders of magnitude from burst to 
burst). The width of the shell can be related to the observed 
duration Tcrb of the GRB through the relation Ao - cTcrb ~ 
10'^Ai2cm for a ~ 30 s burst. This connection between ob- 
served burst duration and shell width holds in the internal shock 
model (Rees & Meszai'os 1994; Nakar & Piran 2002) and pho- 
tospheric models (as discussed, for example, in Giannios & 
Spruit 2007) for the prompt emission. The Lorentz factor of the 
flow is less constrained. From theoretical arguments related to 
the "compactness problem", it is expected to be yo > 100. We 
adopt yo - IQ^'^yi.s as reference value. Among the shell pa- 
rameters, (Tq is the most model dependent quantity and we deal 
with it as a free parameter. In fireball models cro <s; 1 while in 
PDFscro> 1. 

Before studying the initial interaction of arbitrarily magne- 
tized ejecta with the external medium, we summarize the re- 
sults of the unmagnetized cro = case. These results facilitate 
the presentation of the more general case where cro is arbitrary. 

2.2.1. Non-magnetized ejecta 

The deceleration of non-magnetized ejecta has been well stud- 
ied with both analytical (Sari & Piran 1995; Meszaros & 
Rees 1997a) and numerical approaches with 1-dimensional 
(Kobayashi et al. 1999) and 2-dimensional (Granot et al. 2001; 
Cannizzo et al. 2004; Meliani et al. 2007) hydrodynamical sim- 
ulations. The 2-dimensional studies are important to follow the 
late-time lateral spreading of collimated ejecta. Here, we are 
interested in the initial phases of the deceleration of the ejecta. 
These can be followed assuming spherical symmetry since the 
initial Lorentz factor of the flow is likely larger than the in- 
verse of the jet opening angle and, therefore, lateral spreading 
is unimportant. 

The hydrodynamical evolution of the shell-external 
medium configuration is rather well understood. The interac- 
tion of two media leads to a forward shock that propagates into 
the external medium. At the same time a second shock (the 
so called reverse shock) propagates into the ejecta. The for- 
mation of the reverse shock always takes place in cold, non- 
magnetized ejecta since the shocked material moves superson- 
ically with respect to the unshocked shell. This can be for- 
mally shown by solving the initial Riemann problem at the 
ejecta-external medium interface (see, for example RezzoUa & 
Zanotti 2001). We will see that this is not always the case in 
strongly magnetized ejecta where MHD waves can transfer in- 
formation at relativistic speeds. The shell can, thus, smoothly 
adjust its properties upon interacting with the external medium. 

The strength of the reverse shock depends on the ratio of 
the densities of the shell and the external medium and on the 
bulk Lorentz factor of the flow. It can be shown (Sari & Piran 



1995) that the strength of the reverse shock can be conveniently 
parameterized by 
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Aoy^/3' 



(3) 



where I = {3E /Annf^nipC^y^^ is the Sedov length. In the limit 
where ^ » 1 the reverse shock is Newtonian and, if ^ 1, 
relativistic (Sari & Piran 1995). 

The ^ quantity can be related to the more familiar "decel- 
eration radius" rdec and the "spreading radius" r,. The first is 
defined as the radius where the ejecta accumulate from the ex- 
ternal medium a mass y^' times their own mass, resulting in 

2/3 

'"dec = I/Jq ■ The is the radius where a shell with mildly rel- 
ativistic proper motions starts spreading radially = ^qJq (see 
Sect. 3.1). Using the last expressions, ^ can be expressed as 



(4) 



In terms of rdec and rj, the condition for a relativistic reverse 
shock is that rjec /"s ■ 

Using as reference value for the external medium number 
density = 10"«o cm""* we have that 



.1/6 



^ = 0.73 



n}jw,iy,'r 



(5) 



'0 ^12 '2.5 

showing that ^ possesses a very weak dependence with £53 and 
no- For various GRB parameters ^ can typically vary in the 
0.1 $ ^ $ 10 range allowing for both Newtonian and relativistic 
reverse shocks in different bursts. 

In the Newtonian reverse shock limit, the ejecta do not de- 
celerate much during the reverse shock crossing. They enter 
the deceleration phase once they accumulate from the external 
medium mass that is a factor y^ smaller than their own mass. In 
the relativistic reverse shock case substantial deceleration takes 
place already when the reverse shock crosses the ejecta. In both 
relativistic and Newtonian reverse shock limits the radius of the 
reverse shock crossing r^H is given by the expressiorQ (Sari & 
Piran 1995) 



.3 ^l/4 



(6) 



stands for hydrodynamical (non- 



('•s'-dec)' 

where the index "H" 
magnetized) case. 

In the Newtonian reverse shock case, the radius of the re- 
verse shock crossing is modified if one allows for mildly rel- 
ativistic relative speeds within the shell. The relative motions 
lead to shell spreading and strengthening of the reverse shock. 
This effect is considered in Section 3.1 

After the reverse shock reaches the back part of the ejecta, 
there is an initial phase of interaction where shocks and rar- 
efaction waves cross the ejecta. Gradually most of the energy 
is passed in the shocked external medium and the whole struc- 
ture relaxes to the self-similar structure described in Blandford 
& McKee (1976). From this point on, the evolution of the blast- 
wave depends only on the energy E and the density of the ex- 
ternal medium «„. 



' The reverse crossing radius differs in the two limits by a small 
factor which we ignore here. 
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2.2.2. Arbitrarily magnetized ejecta 

The dynamics of the deceleration of magnetized ejecta has 
not been studied in the same detail. To understand the various 
phases of the deceleration of the ejecta and the transfer of their 
energy into the shocked external medium one needs to follow 
the dynamics of ejecta-external medium interaction in detail. 
The internal evolution of the shell is particularly important in 
this respect. The study of the complete deceleration dynamics 
is the goal of the relativistic MHD simulations to be reported 
elsewhere (Mimica et al. 2007). Here we are concerned with 
the very initial phase of the interaction where the reverse shock 
is (or is not) crossing the ejecta and derive the conditions for 
such shock to form. 

The initial phase of reverse shock crossing has been stud- 
ied by Zhang & Kobayashi (2005) by solving the ideal MHD 
shock conditions for arbitrarily magnetized ejecta with a domi- 
nant toroidal field (see also Fan et al. 2004 for the case of mildly 
magnetized ejecta). Their analysis describes the reverse shock 
crossing phase provided that there is a reverse shock fanning. 
They show that the reverse shock crosses faster ejecta of high 
magnetization ctq with respect to un-magnetized ejecta of the 
same E, Aq, Fq. They show that the reverse shock crossing ra- 
dius scales as oc 1 / -sjd\) for ctq s> 1 and can be written 

r,,^P'X^'C,^(r4J'''C,. (7) 

The coefficient Ca depends mainly on the magnetization ctq 
and very weakly on the relative Lorentz factor of the shocked- 
unshocked ejecta (see Fig. 3 in ZK05). For the purposes of this 
work, we adopt the expression Ca ~ 1 / VI + ctq which is rather 
accurate in reproducing their results. Note that this expression 
for Ca gives also the correct zero magnetization limit 

3. Condition for existence of a reverse shock 

Cold, non-magnetized ejecta are always crossed by a reverse 
shock upon interacting with the external medium. This is the 
case since the sound speed of the ejecta is low and does not 
allow for fast transfer of the information of the interaction with 
the external medium throughout their volume. On the other 
hand, in a flow that is strongly magnetized and sub-fast mag- 
netosonic (as in the Lyutikov & Blandford 2003 model) there 
is no reverse shock forming. The flow adjusts gradually to the 
changes of the pressure in the contact discontinuity that sepa- 
rates the magnetized flow from the shocked external medium. 
Here we generalize to arbitrarily magnetized ejecta and derive 
the condition for the formation of a reverse shock. 

An important radius to be considered for the existence of 
a reverse shock is the "contact" radius r^ at which fast MHD 
waves cross the width of the shell. At this radius, the infor- 
mation on the conditions of the front part has been transfered 
to the rest of the shell. After this radius is crossed the magne- 
tized shell comes in pressure balance and adjusts its structure 
on a timescale shorter than the expansion time. Any gradual 
changes of the physical properties at the contact discontinu- 
ity that separates the shell from the shocked external medium 
lead to smooth evolution of the shell properties. The very ex- 
istence of the contact radius implies that difi'erent parts of the 



flow have dropped out of MHD contact during the acceleration 
phase. Our analysis applies for super-fast magnetosonic flows 
(i.e. flows with yo > VI + o"o) for which this is the case. 

A second important radius for the existence of a reverse 
shock is (see eq. (7)). If rc > rrs, a reverse shock forms and 
crosses the ejecta before MHD contact within the shell is es- 
tablished. In this case, internal evolution does not take place in 
the shell before the reverse shock crossing. On the other hand, 
if Tc < Tis, there is no reverse shock forming and does not 
have a physical meaning. The condition for the formation of a 
reverse shock is, thus, to set rc > r^. Note that this condition 
states the obvious fact that there is a reverse shock forming as 
long as it propagates supersonically into the ejecta. As we dis- 
cuss in Section 3.3, a different condition for the existence of a 
reverse shock adopted by ZK05 leads to incorrect estimates on 
the magnetization needed to suppress the formation of a reverse 
shock. 

The "contact" radius is given by the distance the shell 
travels before a fast-magnetosonic wave crosses the shell. 
Any news from the front travels into the shell with speed 
/3a = Vo"o/(1 + ctq) (with a corresponding Lorentz factor 

- VI + o"o) in the frame comoving with the shell. As seen 
in the central engine frame, this wave moves with Lorentz fac- 
tor y„, given by the expression ja ^ (yw/yo + yo/TH)/2. The 
last expression is accurate for jo, y,,. » 1 and can be solved for 

resulting in -/„ = yd VI + o-q - -\fd\)). The shell comes into 
contact when the wave crosses it. This takes place at radius 




where we have used that vv/c ^ 1-1/2-)/^, and vq/c ^ l-l/2y^. 

The contact radius depends only on the shell properties and 
not on those of the external medium. The reason is that the 
speed of propagation of fast-magnetosonic waves depends only 
on the magnetization of the shell as a function of the radius cr(r) 
which turns out to be equal to the initial value ctq. To check 
such an assertion, one needs to consider that a shell that ex- 
pands spherically decreases both its density and its magnetic 
energy density as r^^. Therefore, cr(r) remains uniform in the 
shell and equal to its initial value (eq. 1) until, at least, the shell 
reaches r^. 

From eq. ^ we see that in the limit of cro » 1, rc ^ 
AoyQ/2cro oc cTq' while for ctq — > 0, r^. — > oo. In the hydro- 
dynamical (ctq = 0) limit rc > ris and a reverse shock crosses 
the ejecta at rrs. On the other hand, since the contact radius 
scales with the magnetization as rc oc cTq', it is reduced faster 

— 1/2 

than the rjs oc cr^ for increasing ctq. There is always a mag- 
netization for which rrs and rc become equal and the reverse 
shock is suppressed. 

The condition for the existence of the reverse shock is given 
by setting r^ < rc and using eqs. (|7]i, dHJ and the definition ^ 
of ^. The line that separates the reverse shock from the non- 
reverse shock regime in the ^-ctq parameter space is 

= ^3 ^ ^ 2cr„ - 2 Vcrod+cro)). (9) 

\ rs / cro 
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Fig. 1. Regime of existence of a reverse shock in the ^-ctq plane. 
The ^ is defined by eq. (O (or eq. (4)) and ctq stands for the 
Poynting-to-kinetic flux ratio of the flow at the onset of the af- 
terglow. Ignoring radial spreading of the ejecta, the lower left 
part with respect to the black, dashed line shows the regime 
where there is a reverse shock forming for a constant-density 
external medium. For ctq > 1 there is no reverse shock form- 
ing for a large parameter space that describes the properties of 
the GRB ejecta. The solid line shows the frontier between "no 
reverse shock" and "reverse shock" when radial spreading of 
the shell is taken into account. In the case of shell spreading, 
the shaded area is characterized by a mildly relativistic reverse 
shock (see Sect. 13. lb . Similarly, the red lines show the frontier 
between "no reverse shock" and "reverse shock" in the wind- 
like external medium. 

In the limit of small and large cTo the last expression simplifies 
respectively to 

^ = -jj^ for a-Q <s; 1, 



(10) 



1 



(4cro)i/3 



for (To » 1. 



Although, the expression ^ = (4cr())"'/^ is formally derived 
for (To » 1, it is still a fair approximation for (To ^ 0.1 and will 
be used in the following discussion. Solving the last expression 
for (To, and using the reference values of the GRB parameters 
we find that a reverse shock forms when 



(11) 



The curve defined by eq. (|9]l is shown in Fig. 1 with a 
dashed line that turns into solid for (To ^ 0.3. It separates the ^- 
cTy parameter space into regions with the upper right one char- 
acterized by no reverse shock. In the (To <K 1 regime the ^ - 1 
dash-dotted line separates the Newtonian reverse shock (^ > 1) 
regime from the ^ < 1 relativistic reverse shock one. Note that 
a small magnetization suffices to suppress a reverse shock for 
large ^ and vice versa. 

As it is obvious from eq. (fTTT i. for typical GRB parameters 
for a strongly magnetized flow with ctq > 1 there is no reverse 



shock forming. In this case, the onset of the afterglow emission 
will not be characterized by an optical flash (Meszaros & Rees 
1997a; Sari & Piran 1999) that is expected to be the signature of 
the reverse shock emission. Note, however that this conclusion 
depends sensitively on some of the GRB parameters. This is 
especially true for the Lorentz factor of the GRB flow which 
is also, perhaps, the least constrained parameter Had we used 
a reference value of yo - 1000 for it then the transition to no 
reverse shock would take place for (To ~ a few tens. This means 
that one can still have a reverse shock forming in a yo ~ 1000, 
(To » 1 flow. However, when a reverse shock forms in a high 
(To flow it is not expected to result in bright emission (at least in 
the ideal MHD limit discussed here) since the shock dissipates 
only part of the kinetic energy which is a small fraction of the 
total energy of the shell (ZK05). 

So far, we have focused in ejecta with initially uniform mo- 
tion, i.e. with no radial spreading. Shell spreading introduces 
new features to the picture presented so far These features are 
the topic of the next section. 

3.1. Shell-spreading effects 

It is quite likely that the ejecta are characterized by relative mo- 
tions after the acceleration, collimation and internal interaction 
phases. These relative motions may well be mildly relativistic 
in the shell rest frame leading to ejecta with moderate varia- 
tions of the Lorentz factor 6yo ~ yo along their width. 

Ejecta with their front part moving faster than their rear 
one with moderate variation of the bulk Lorentz factor start 
spreading at radius r, ^ Tq^o (Sari & Piran 1995; Kobayashi 
et al. 1999). Because of the spreading, the width of the shell in- 
creases linearly with distance at larger radii. It has to be noted 
that Ao/c does not necessary correspond to the time during 
which the central engine is actively releasing energy into the 
outflow. As pointed out by Aloy, Janka & Miiller (2005), in 
case of short GRBs, the initial shell width A^ - cTg ija being 
the time during which the central engine actively releases en- 
ergy) may stretch radially by a factor ~ 2 - 20 due to the same 
mechanism invoked for shell spreading at later times, namely, 
the relative velocity difference between the forward and rear 
radial edges of the ejecta. 

In the case of non-magnetized ejecta with ^ > 1, the shell 
starts spreading before the reverse shock crossing. Radial ex- 
pansion leads to density decrease in the shell which is faster 
than 1/r^. As a result, by the end of the reverse shock cross- 
ing the shock is mildly relativistic. When ^ < 1, the reverse 
shock crosses the shell before it reaches the spreading radius rs- 
The relativistic reverse shock case is, hence, not affected by the 
spreading of the shell since, by the time the shock has crossed 
the shell the density of the latter has not dropped appreciably 
in comparison to the non-spreading case (Sari & Piran 1995). 
One may also express the condition of formation of a relativis- 
tic reverse shock in terms of the deceleration radius rdec and r^, 
which reads rdec ^ /"s- 

Spreading of the shell influences the interaction of the 
ejecta with the external medium for magnetized ejecta as well. 
If the contact radius is larger than the spreading radius, the 
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shell spreads before fast MHD waves produced in the front 
part of the shell catch up with the rear part. Using the esti- 
mate Vs - TgAo for the spreading radius and eq. (|8j for r^, one 
finds that < for cry < 0.3. For low enough ctq the shell does 
not come into MHD contact because of radial spreading. Under 
these conditions, a mildly relativistic RS forms just like in the 
ctq -Q case. 

The shaded region of Fig. 1 shows the regime where spread- 
ing of the shell leads to the formation of a reverse shock. Quite 
generally for erg < 0.3, there is always a reverse shock form- 
ing in a shell that spreads because of mildly relativistic proper 
motions. 

3.2. Deceleration in a wind environment 

A possible environment of the deceleration of the GRB flow 
is that of stellar wind (i.e. "collapsar" model for long GRBs; 
Woosley 1993). The density profile of the external medium in 
this case is - Ajr^. Here, we derive briefly the condition for 
the existence of a reverse shock in such density profile. 

The reverse-shock crossing radius for un-magnetized ejecta 
that decelerate in wind density profile is rj^^,, ~ (rs'"dec,w)''^ 
(Sari & Piran 1995), where the deceleration radius is defined as 
''decw - Ivi/yf) - E / A-nAmpj^c^ . The is the Sedov length for 
the wind density profile. If the ejecta is magnetized, the reverse 
shock crosses them faster with respect to the hydrodynamical 
case. In the high-cro limit the (reverse) shock conditions lead 
to rrs,w - '"re w/"^!) (^'g- ZK05). An approximate expression for 
the radius of the crossing of the reverse shock that reproduces 
both the ctq — > and ctq » 1 limits is 



l+cr,, 

The condition for existence of a reverse shock is that the rre.w < 
rc. The contact radius depends only on the shell properties and 
not on the density profile of the external medium; it is defined 
by eq. ([S). The line that separates the "no reverse shock" from 
the "reverse shock" regime is found by setting rre,w = rc. It 
results in: 

The line defined by eq. (13) is shown in the ^-ctq plane of Fig. 1. 

If the shell is characterized by mildly relativistic proper mo- 
tions, its spreading radius is rs ~ Aoyg. Just as in the case of 
constant density external medium, spreading of the shell leads 
to the formation of a reverse shock when ctq < 0.3. 

3.3. Non-ideal MHD effects 

In this work, we focus on the existence of a reverse shock 
(and the associated emission) coming from the interaction of 
a shell, permitted by a large scale magnetic field, with the ex- 
ternal medium. Furthermore, this analysis is based on the ideal 
MHD shock conditions derived by ZK05. However, non-ideal 
MHD efi'ects may also be of relevance in the flow. 



It is possible, for example, that the magnetized flow is 
launched with magnetic fields that contain a random com- 
ponent; one that changes polarity on short scales in the ra- 
dial direction (Thompson 1994; Drenkhahn & Spruit 2002; 
Thompson 2006). Furthermore small scale field reversals may 
develop during the acceleration phases as result of MHD insta- 
bilities in the jet (Giannios & Spruit 2006). Unless magnetic re- 
connection is efficient in dissipating the random component of 
the magnetic field (e.g. Drenkhahn & Spruit 2002), small scale 
reversals of the field along the radial direction are still present 
in the afterglow regime. Such field structure can have an ef- 
fect on the shock conditions and therefore on the condition for 
the existence of a reverse shock (for discussion see Thompson 
2006). 

Dissipation of part of the magnetic energy that is contained 
in the flow at large radii can result to particle acceleration and 
emission with a variety of observational consequences. This 
emission is not directly related to shock-accelerated particles. 
The reverse shock can, for example, lead to forced reconnection 
of a reversing magnetic field. The released energy can power 
the prompt GRB emission (Thompson 2006). Magnetic dissi- 
pation can also take place during a later stage of the deceler- 
ation of the ejecta. In this picture, slowing down of the ejecta 
leads to the revival of MHD instabilities in regions of the flow 
that come into causal contact again. Such delayed dissipation 
can power the afterglow flares (Giannios 2006a) observed in 
many early X-ray afterglows (Burrows et al. 2005; Zhang et al. 
2006). 

3.4. Comparison witli previous work 

The early phases of interaction of the GRB ejecta have been 
previously studied in the hydrodynamical limit (i.e. erg — > 0; 
Sari & Piran 1995) and for a sub-fast magnetosonic, strongly 
magnetized, flow (i.e. Tg < 1 +crg; Lyutikov & Blandford 2003; 
Genet et al. 2006). The outcome of those studies is that in the 
former limit there is always a reverse shock developing, in the 
latter the flow adjusts smoothly to the properties of the shocked 
external medium and there is no reverse shock. Those results 
are in agreement with these reported here. 

The intermediate range of magnetization has been studied 
in ZK05 by solving for the shock conditions for the forward 
and the reverse shock. Their study is applicable for the param- 
eter space where there is a reverse shock forming. They also 
estimate the conditions for the formation of a reverse shock 
assuming a constant density external medium. They find that 
large values of <jq ~ 100 are needed (for the reference GRB pa- 
rameters) to suppress the existence of a reverse shock. This es- 
timate is more than two orders of magnitude different from the 
one presented here. To probe the difference, we briefly sketch 
the ZK05 condition for the existence of a reverse shock. 

ZK05 define the deceleration radius as the radius where 
the ejecta accumulates from the external medium mass a fac- 
tor 1/yg less than the external medium mass: r^^ - l/y^^^il + 
(To)'^-'. The conditions for the existence of a reverse shock that 
they use are that (i) the flow speed exceeds the fast magne- 
tosonic speed (i.e. yo > VI + o"o) and (ii) the magnetic pres- 
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sure of the shell is lower than that of the shocked external 
medium at the deceleration radius r^^: i.e. that B^J^ny^^ < 
(4/3)7^eC^. In the last expression, one uses that the comov- 
ing magnetic field is Bq.co = ^o/To ™d the forward relativistic 
shock condition that gives for the gas pressure of the shocked 
external medium: Pg = (4/3)y^eC^- Using eq. (2) and fo- 
cusing in the ctq » 1 limit, the condition for ctq for the ex- 
istence of a reverse shock can be written ctq < 190/^^. The 
boundary that separates the shock from the non-shock regime 
is (TzK - 190/^^. It is very different with respect to that of 
eq. ([Tot derived here. 

What is the source of this discrepancy? ZK05 apply their 
condition in the deceleration radius and not the reverse shock 
crossing radius which we believe is the relevant one so as to 
check the existence of a reverse shoclfl Furthermore, the con- 
dition itself for the existence of the reverse shock used in ZK05 
is questionable. For magnetization ctq ~ ctzk the "contact" ra- 
dius Tc is much shorter than the "deceleration" radius r^^ and 
the reverse shock crossing radius ris. This means that one gets 
the erroneous result that the reverse shock moves subsonically 
with respect to the unshocked material. The fact that < rf^'^ 
indicates instead that much before the flow reaches r^^ the 
shell has already evolved. The spherical geometry starts to play 
an important role and the problem becomes that of dynamical 
evolution of the shell properties and not just shock conditions 
of a non evolving shell. 

On the other hand, when a reverse shock does form then 
it crosses the shell before its different parts have time to come 
into contact through exchange of MHD waves. This means that 
the properties of the shell upstream to the shock do not evolve 
much before the shock crossing. In this case the analysis of 
ZK05 is applicable. This justifies the use of their result on the 
reverse shock crossing radius given by eq. (|7|i- 

4. Discussion 

The role of magnetic fields in GRB flows remains uncertain. 
In fireball models, the GRB flow is launched with dynamically 
unimportant magnetic fields while MHD energy extraction can 
lead to a Poynting-flux dominated flow. Observations of the 
prompt emission alone appear hard to distinguish between fire- 
ball or PDF modelfl Early afterglow observations can, how- 
ever, probe the magnetic content of the ejecta. The magnetiza- 
tion can leave its signature in the early forward shock (Lyutikov 
& Blandford 2003; Genet et al. 2006) and the reverse shock 
emission (Kumar & Panaitescu 2003; Zhang et al. 2003; Fan et 
al. 2004; ZK05). 

4.1. Magnetization and existence of a reverse shock 

Here, we derive the condition for the existence of a reverse 
shock in the ejecta that contain a large scale magnetic field, 

^ Note also that the definition of the deceleration radius is contro- 
versial (see Lyutikov's 2005 comments) 

^ with the possible exception of multi GeV prompt emission that 
may show signature of pair cascades in a neutron rich fireball that is 
absent in PDF flows; Koers & Giannios (2007). 



assuming the ideal MHD description is adequate. The ejecta 
decelerate interacting with the circumburst medium. External 
medium with both constant and wind-like density profiles are 
considered. It is shown that, for high enough magnetization, the 
ejecta interact smoothly with the external medium and deceler- 
ate without the formation of a reverse shock. The conditions for 
the formation of a reverse shock are derived by demanding that 
the reverse shock propagates supersonically in the ejecta. The 
conditions are summarized by eqs. ^ and (fTTT i and Fig. 1. Our 
results are substantially diff'erent from those of ZK05 who do 
not consider the internal evolution of the ejecta when deriving 
their criterion for the existence of a reverse shock. 

One important finding of the study is that ctq > 1 suffices to 
suppress the formation of a reverse shock for a large parame- 
ter space that characterizes GRB flows. This parameter space 
is larger for a constant density external medium in comparison 
to the wind case. This means that only a small fraction of PDFs 
are expected to show reverse shock emission. Even when a re- 
verse shock forms in a high ctq flow, it is not expected to result 
in bright emission (at least in the ideal MHD limit discussed 
here) since the shock dissipates only part of the kinetic energy 
which is a small fraction of the total energy of the shell (see, 
e.g., ZK05). On the other hand, allowing for radial spreading of 
the ejecta, practically all the flows with ctq < 0.3 are expected to 
be crossed by a mildly-relativistic or relativistic reverse shock. 

This work contributes in understanding the early deceler- 
ation dynamics of magnetized ejecta. The complete evolution 
of the system is not considered here. Important questions re- 
lated to the timescale over which the magnetic energy of the 
shell (which may dominate the total energy budget) is passed 
into the shocked external medium remain unexplored. We are 
currently attacking these questions through relativistic MHD 
simulations to be reported elsewhere (Mimica et al. 2007). 

4.2. Connection to observations 

A clear prediction of the fireball models is that there is a re- 
verse shock crossing the ejecta on the onset of the interac- 
tion of the ejecta with the external medium. The reverse shock 
emission is expected to result in a short lived, bright emis- 
sion event. Although optical flashes have been observed in a 
few bursts (e.g. GRB990123, Akerlof et al. 1999; GRB02121 1, 
Fox et al. 2003) and are believed to be powered by the reverse 
shock emission, they are rare events. This may be a result of a 
more luminous forward shock masking the reverse shock opti- 
cal emission (McMahon et al. 2006; Mundell et al. 2007b). 

Our analysis points to a different possibility, i.e. that the 
paucity of bursts with clear reverse shock emission is a result 
of the frequent reverse shock suppression. It is shown here that 
this happens for a large GRB parameter space in a strongly 
magnetized ffow. Moreover, in cases where reverse shock emis- 
sion is identified, modeling of the microphysical parameter 
connected to the magnetic field strength indicates ejecta with 
dynamically important magnetic fields (Fan et al. 2002; Kumar 
& Panaitescu 2003; Zhang et al. 2003; ZK05) in agreement 
with our suggestion. 
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Polarized emission, result of synchrotron emitting, (re- 
verse) shock-accelerated particles is expected from a coherent 
magnetic field. Note, however, that the lack of significant op- 
tical polarization in the early afterglow does not necessarily 
show the absence of large scale field in the ejecta but may be 
result of the absence of reverse shock (suppressed by the strong 
field). This is particularity likely if there is no evidence for the 
distinct emission patterns expected from a reverse shock. In 
light of these considerations the claim of Mundell et al. (2007a) 
that the low early polarization signal in GRB060418 rules out 
the presence of large scale magnetic field in the ejecta is pre- 
mature. 

This study has focused in the ideal MHD limit. On the 
other hand non-ideal MHD eff'ects can lead to additional energy 
release and emission during the afterglow phases of strongly 
magnetized flows. This emission is not directly related to shock 
accelerated particles and therefore to the existence or not of a 
reverse shock. On the other hand, the crossing of the reverse 
shock can lead to forced reconnection in the ejecta that powers 
the GRB emission (Thompson 2006). Furthermore, the decel- 
eration of the ejecta can lead to revival of MHD instabilities 
that lead to magnetic dissipation. Such late-time energy release 
can power the afterglow flares observed in many early X-ray 
afterglows (Giannios 2006a). 
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